A unique parametrization of the shapes of secondary dilepton 
spectra observed in central heavy- ion collisions at CERN-SPS 

energies 



K. Gallmeister", B. Kampfer", O.P. Pavlenko ' 6 , C. Gale c 

a Forschungszentrum Rossendorf, PF 510119, 01314 Dresden, Germany 

^Institute for Theoretical Physics, 252143 Kiev - 143, Ukraine 

^ [ Thysics Department, McGill University, Montreal, QC, H3A 2T8, Canada 

> 
O 

00 '. Abstract 

^ ! A unique parametrization of secondary (thermal) dilepton yields in heavy-ion ex- 
periments at CERN-SPS is proposed. This parametrization resembles a thermal qq 
annihilation rate. This is inspired by the observation that lepton pair production rates 

O ■ are quantitatively similar, whether expressed in a hadronic or partonic basis. Adding 



the thermal yield and the background contributions (hadronic cocktail, Drell-Yan, cor- 
related semileptonic decays of open charm) the spectral shapes of the CERES/NA45, 
NA38, NA50 and HELIOS/3 data from experiments with lead and sulfur beams can 
be well described. 
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1 Introduction 

The typical temperature scales in heavy-ion collisions at CERN-SPS energies, extracted 
from hadron abundances M and hadron transverse momentum spectra 0, [| |J, are in 
the order of T c (or somewhat less), where T c stands for the expected deconfmement 
and/or chiral symmetry restoration temperature in strongly interacting matter. For 
temperatures around that value, a critical comparison of dilepton-producing channels 
suggests that, in the intermediate invariant mass region, the sum of various hadronic 
processes essentially equals the quark-antiquark annihilation rate into lepton pairs 
||. Similarly, at lower invariant masses the various contributions to the lepton pair 
spectrum add up to produce a structure where the p meson has acquired a considerable 
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width ||. Remarkably, there also the net dilepton spectrum closely mimics that from a 
free gas of annihilating quarks and antiquarks 0. It is tempting to ascribe this duality 
of representations to a full or partial restoration of chiral symmetry. This interesting 
conjecture however needs yet to be put on a firm theoretical basis. We simply intent 
here to exploit its empirical power. The interested reader can consult references || [| 
for early analyzes, for a recent review, and [|T0| , |TT| , |T2]j for a recent application. 

We employ here the qq rate as a convenient parametrization of the dilepton emissiv- 
ity and analyze uniquely both the e~e + and channels of the dilepton spectra, 

obtained in sulfur and lead beam reactions at the CERN-SPS. We shall make use of 
the analytical simplification brought about by only considering this simple process, to 
attempt a global fit to the world data of lepton pair production in central heavy-ion 
collisions at SPS energies. Since we use the qq rate in the full invariant mass range 
considered and for the full time evolution of matter as well, one could phrase our ap- 
proach as a test of an "extended duality" hypothesis. Again, for the moment we refrain 
from attempting to formulate a theoretical foundation for this fact and we first seek to 
check its validity in a study of heavy-ion phenomenology. 

The dilepton spectra are obtained by a convolution of a rate with the space-time 
evolution of matter. The time evolution is split up in several stages which can be dealt 
with separately: (i) First, on very short time scales there are hard initial processes 
among the partons, being distributed according to primary nuclear parton distribu- 
tions, such as the Drell-Yan process (to leading order qq — > 11) and charm production 
(to leading order mainly gg — > cc). (ii) On intermediate times scales there are the 
so-called secondary interactions among the constituents of the hot and dense, strongly 
interacting matter. This stage is often denoted as thermal era and the emitted dilep- 
tons as thermal dileptons. (iii) If the interactions among the hadrons in a late stage 
cease, there are hadronic decays into dileptons and other decay products. 

In the present work we focus on a parametrization of the dilepton yield from stage 
(ii) with a minimal parameter set and assume that the background contributions from 
stages (i) and (iii) are known, either by up-scaling the Drell-Yan and open charm 
yields from pp to AB collisions or by using directly the experimentally determined 
hadronic cocktail. One should be aware that such a schematic description may suffer 
from several deficits, such as missing secondary Drell-Yan processes ||13|| , or hadronic 



final state interaction of open charm |L4| or non-equilibrium effects ||15|| . Also, if the 
temperature is initially large enough, i.e. T ^> T c , and the matter resides in a deconfined 
state, one expects a dilepton spectrum being different from the simple qq rate both in 
the low-mass region [|16j and in the intermediate- mass region [17] (cf. Hl8|). 
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The aim of the present work is to look for a minimal description of data in a similar 
spirit as one comfortably parametrizes the transverse momentum spectra of hadrons 
by effective slope parameters and afterwards maps on flow and freeze-out temperature 
and non-equilibrium components. One can also draw from experiences with direct 
photons, where a simple parametrization of hadronic rates has simplified calculations 
considerably |[T9| . 

Our paper is organized as follows. In section 2 the model for the dilepton emission 
rate is presented. The analyses of the experimental spectra are performed in section 
3. Flow effects and time evolution effects are considered in section 4. The results are 
summarized in section 5. The calculation of the background contributions is explained 
in the appendix A. 



2 Model parametrization 

Let dN/(d 4 x d 4 Q) be the Lorentz invariant dilepton production rate from matter in 
local thermal equilibrium being characterized by a temperature T, chemical potentials 
fii, and four- velocity u v of the flow. We base the production rate 0, on the lowest- 
order quark- ant iquark annihilation rate 

dN 5a 2 ( u-Q) 

exp —}, (1) 



d 4 Q d 4 x 36vr 4 * I T 

where = (M± coshY, M± sinhY, Q±) is the lepton pair four-momentum with 
transverse mass M±, transverse momentum Q±, related to the invariant mass via 



M = y M]_ — Q 2 ±1 and rapidity Y . The above rate is in Boltzmann approximation, 
and a term related to the chemical potential(s) is suppressed. The space-time inte- 
gration can be performed if a model for the set of parameters T, /ij, and their 
space-time dependence is at disposal. An obviously strong approximation is to replace 
all state variables by averages 

^-^expj-^j [d 4 x (2) 

It is one of the purposes of this paper to investigate the validity or the consequences 
of this potentially extreme, but remarkably simple assumption. In section 4 we discuss 
that the flow effects can be neglected in some region of the phase space. Then, by 
assuming a thermal source at midrapidity F C ms, one can make the replacement — > 
(cosh y cms , sinh y cms , 0). Taking the space-time volume / d A x = N eS as normalization 
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factor (which now can also include effects of finite chemical potetials) one gets the 
approximation 

dN 5a 2 AT f M ± cosh(F - y cms ) \ 

where T eff = (T). The two parameters T eS and iV eff are to be adjusted to the experi- 
mental data. 

In a more detailed description these two parameters are mapped on a much larger 
parameter space, which however is constrained by hadronic observables and allows a 
detailed microscopic justification of T eS and N e g. 

In what follows we use Eq. (3) and confront it with the data. In doing so we implement 
the corresponding detector acceptance, which is most conveniently done by generating 
the six-fold differential rate 

dN 1 dN 

= , (4) 

P±idp±ip±2dp± 2 dy 1 dy 2 d(pid(p2 2n d 4 Q } 

where p_i_i,2, 2/1,2 and 4>i,2 denote the transverse momenta, rapidities and azimuthal 
angles of the individual leptons 1 and 2, which must be appropriately combined to 
construct the pair mass M, the pair transverse momentum Q± and transverse mass 
M, . 



3 Analysis of dilepton spectra 
3.1 Lead beam data 

The dilepton experiments in the e~e + and channels with the lead beam at CERN- 
SPS use as targets either Au or Pb. We neglect the small differences of the target nuclei 
and attempt a unique parametrization. The rapidity coverage of the two experiments 
is also fairly symmetric around midrapidity for these symmetric collisions. 



3.1.1 CERES experiment 



A comparison of the model defined by Eqs. (3, 4) with the lead beam data [2C] is 
displayed in Fig. [l|. With the parameter set of F cms = 2.9, T eff = 170 MeV and N cS = 
3.3 x 10 4 fm 4 a fairly good description of the data is accomplished. We mention that the 
use of external information, namely the hadronic decay cocktail and the normalization 
of (N c h) = 250, is essential for describing the CERES e + e~ data |20| in central reactions 
Pb(158 AGeV) + Au. The normalization to (N c h) = 250 corresponds to a centrality 
criterion of cr triggeT / a tot = 0.3 pO . 
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We model the CERES acceptance by p ±1 , 2 > 0.2 GeV, rft$ = 2.1 •• -2.65 (with rj 
as pseudo-rapidity) and a relative angle 0^ b > 35 mrad between the electron and 
positron. 



In agreement with other models [21, 22] we also describe the Q± dependence in various 
M bins (see Fig. |). 

It should be emphasized that the above parameters T eS and N eS deliver also an 



optimum description of the direct photon data 123], as shown in 24 



3.1.2 NA50 experiment 



For a description of the NA50 /x 4 ' \i~ data [25] in central reactions Pb(158 AGeV) + 
Pb, the Drell-Yan contribution and the correlated semileptonic decays of open charm 
mesons are needed. The latter ones are generated with PYTHIA |26| where the Drell- 
Yan K factor of 1.5 is adjusted to the data [27, |28[] (cf. Appendix |A.1|) and the open 
charm K factor of 5.7 to the compilation cross sections of identified hadronic charm 
channels [2S] and [i + \i~ data in the reaction p(450 GeV) + W |3(J (cf. Appendix |A.2|) . 

To translate the cross sections delivered by PYTHIA into rates we use a thickness 
function of 31 mb _1 for central collisions Pb + Pb. Actually, however, the centrality 
criterion is a selection of data from E? bin 9 with impact parameter average (b) < 3.3 
fm and a participant number (N paxt ) = 381 ± 7 according to [p5[] . 

The NA50 acceptance is modeled by Y lab = 2.9 • • ■ 3.9, a Collins-Soper angle 
| cos 0cs| < 0.5 and the condition of E > E m i n for the minimum energy of a muon in the 
laboratory, where E min = E + AE with E = 11.5 GeV and AE = 16000(6 - 0.065) 2 
GeV, 0, 13000(0 - 0.090) 2 GeV for 0.037 < 6 < 0.065, 0.065 < 6 < 0.090, 0.090 
< 6 < 0.108, respectively. Note that these cuts describe the NA50 acceptance only 
approximately. 

The resulting invariant mass and transverse momentum spectra, including the ther- 
mal source contribution, are displayed in Fig. |2|. The thermal source, with parameters 
T fj and N e s adjusted to the above CERES data, is needed to achieve the overall 
agreement with data. This unifying interpretation of different measurements has to be 
contrasted with other proposals of explaining the dimuon excess in the intermediate 



mass region either by final state hadron interactions |14| or by an abnormally large 
open charm enhancement P5 |. The latter one should be checked experimentally [|3T 



thus attempting a firm understanding of dilepton sources. Notice that our minimum 
parameter model describes the data equally well as the more detailed dynamical models 

§0. 
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3.2 Sulfur beam data 

Let us now turn to the older sulfur beam data (cf. for a recent survey on models 
and an extensive reference list). Since a much larger rapidity interval is covered (see 
Fig. |3p we smear the source distribution (3) by a Gaussian function with a width of 
a = 0.8, i.e. 

dN 5a\ T f dY' ( M ± cosh(Y-Y') \ j (Y'-Y cms ) 2 \ 

[^eff / r- — = exp { \ exp \ — \ . (5) 



d 4 Q 367T 4 J y/^TTO 2 { T cS j ^ { 2<7 2 

For Y cms we choose 2.45 as suggested by an analysis of the rapidity distribution of 
negatively charged hadrons from NA35 ||3^] as displayed in Fig. |j. The width of the 



dilepton source, a, is somewhat smaller than the width of the hadron distribution, cf. 
Fig. [|. Here we neglect also the small differences of the various target nuclei (Au, U, 
W) and attempt a unique parametrization of the e~e + and channels. 

3.2.1 CERES experiment 

For the CERES e + e~ data |33| in central S(200 AGeV) + Au reactions, corresponding 
to = 125, the published hadronic cocktail is used. The acceptance is the same 

as described in subsection |3.1.1| . The comparison of our calculations with the data is 
displayed in Fig. [5[ A good description of the data is achieved by T e ff = 160 MeV. One 
observes in the left pannel of Fig. [5] indeed a fairly well reproduction of the spectral 
shapes for the choice of the normalization factor N c g = 11.2 x 10 4 fm 4 . (Notice that 
this normalization factor is ununderstandably large. We focus here, however, on the 
shape of the spectra and do not attempt a change of our simple parametrization to 
resolve this issue, e.g. by employing another rapidity distribution. In this context we 
mention the second reference in [^] where, within the same transport code, the lead 



beam data |20| are satisfactorily reproduced but the sulfur beam data |33[] strongly 



underestimated. Similarily, if we use the same normalization factor, as adjusted to the 
HELIOS/3 data (see subsection 3.2.3 below) our resulting spectrum is below the data 
and even outside of the systematical error bars at M ~ 400 MeV, see right pannel of 
Fig. |.) 



We mention additionally that adjusting the normalization to the HELIOS/3 data [37| 



the published upper bounds of the direct photon yields 0] are just compatible with 



our model calculations when adopting the model described in section 2 for photons p4[ 
(see Fig. ||). In contrast, when choosing the normalization, which delivers an optimum 
description of the CERES data, one would be above the upper photon bounds in |34 
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3.2.2 NA38 experiment 

The acceptance of this dimuon experiment is described by Y[ a b = 2.8 •• -4.1 and 



|cos6cs| < 0.5. The comparison of our calculations with the NA38 data for 
the reaction S(200 AGeV) + U in the intermediate-mass and high- mass region is dis- 
played in Fig. [7|. Adjusting the normalization factor to the data we achieve an optimum 
description by N e g = 1.3 x 10 4 fm 4 (left pannel), while the use of the normalization, 
adjusted to HELIOS/3 data, results in a unsatisfactory data description (right pannel 
in Fig. [7]). It should be emphasized that all other data sets we are analyzing are for 
more restrictive central events. It is therefore clear that the required normalization 
factor for the NA38 data is smaller. 

The available transverse momentum spectrum of NA38 is also nicely reproduced in 
shape (see Fig. ^). 

Recently, the NA38/50 collaboration published also dilepton data in the low-mass 
region ||36|| . Since we have no reliable background contribution (hadronic cocktail) 
within the given acceptance at hand we discard the inclusion of these data in our 
analysis. 

3.2.3 HELIOS/3 experiment 

While we can nicely reproduce the Drell-Yan background for the HELIOS/3 experiment 
(cf. |37fl), our PYTHIA simulations deliver another open charm contribution than the 
one used in previous analyses 0. Since the accurate knowledge of the background 
contributions is necessary prerequisite, we use therefore for our analysis the difference 
fi + fi~ spectra of S(200 AGeV) + W and p(200 GeV) + W reactions f37j, thus hoping 



to get rid of the background since these spectra are appropriately normalized. The 
centrality selection for these data is a multiplicity larger than 100 in the pseudo-rapidity 
bin 3.5 • ■ ■ 5.2 resulting in an averaged multiplicity of 134.6 



The acceptance is described by M± > max ^4(7 - 2Y, ab )GeV, y^m^ + (^|^) 2 | 

and the data are binned in the rapidity intervals Yj a b = 3.0 ■ • ■ 3.9, 3.9 • • • 4.4, 4.4 • ■ ■ 7.0. 
It is in particular this wide rapidity span and the binning which require the Gaussian 
smearing of the dilepton source. Considering the rapidity-integrated yield, as done in 
most previous analyses when only the integrated data were at disposal, discards an 
important information and allows a unique description of the CERES lead beam data 
and HELIOS/3 sulfur beam data without any problem |38| . 
As seen in Fig. |9] the difference spectra in the various rapidity bins are fairly well 
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described for N c g = 5.3 x 10 4 fm 4 . (Notice that, in comparison with the lead beam 
data, also this normalization is quite large.) The M±_ spectra in the two available M 
bins for the integrated rapidity bin Y = 3 • • • 4.4 are also described in gross features, 
see Fig. |10|. (Integrating the experimental M± spectra in this figure one finds a factor 
of 2 difference to the spectra in Fig. [] when integrating these over the corresponding 
M intervals. This is accounted for in the Fig. [T0| .) 



4 Discussion of flow and time evolution effects 
4.1 Flow effects 

As qualitatively discussed in [TTJ flow effects affect mainly the Q± or M± spectra. Ex- 



plicit formulae for the dilepton spectra are given in [|TT|] for (i) spherical expansion and 
(ii) longitudinally boost-invariant expansion superimposed on transverse expansion. 
In Fig. [TT] we show the result of Monte Carlo simulations of the rates dN/ dM dY and 
dN/dM± dY at Y = for various values of the transverse flow velocity. One observes 
that the invariant mass spectra, delivered by a Monte Carlo procedure for generating 
the distribution Eq. (4), are indeed insensitive against flow (see Fig. [IT], left panel). 
This is known for case (ii) for some time P5| |. In contrast, the Q± spectra are sensitive 



against flow as demonstrated in the right panel of Fig. |TT|, in particular in the large-Q^ 
region. However, as seen in Figs. [IJ and ^, the large Q± region is dominated by the 
background contributions (hadronic cocktail or Drell-Yan). Therefore, one can indeed 
neglect the flow. Note that the key for this statement is a constraint on transverse 
flow from hadronic data. The analyzes of the m± spectra of several hadron species 
point to flow velocities (v±) in the range from 0.43 c 0, fl) to 0.55 c ||. Since the flow 
is expected to increase continuously with time due to acceleration from pressure, the 
transverse flow of hadrons at kinetic freeze-out is a temporal maximum value. The 
spatial and timewise average of the flow, relevant for dilepton and photon emission, 
must be smaller than the quoted values. 

The situation may change at RHIC and LHC, where higher initial temperatures, and 
consequently larger pressures, cause a stronger transverse flow which could be manifest 
in dilepton Q± spectra. 
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4.2 Time evolution effects 



Let us now discuss time evolution effects. As an example we show in Fig. ^ such initial 
temperatures and the corresponding final temperatures as a function of an normaliza- 
tion factor, N, which deliver the same invariant mass spectra in the range covered by 
the CERES and NA50 lead beam experiments displayed in Figs. [1] and 0. In calculating 
the curves in Fig. |l2]the time evolution according to T = Too + (^Initial — Too) ex P {~ 
V{t) = iV|±gexp {^}, with = 110 MeV, t x = 5 fm/c, t 2 = 8 fm/c is used as in 
[ |TT| 1 (cf. [50]), where A, S are the mass numbers of the colliding nuclei, no = 0.17 fm -3 . 
The temperature evolution starts at Tj n i t i a i an d is stopped at Tfi na i. (Such a tempera- 
ture and volume evolution have been used first in to show that the superposition 
of the thermal qq annihilation rate and the hadronic cocktail reproduces the shape of 



the CERES data obtained with the sulfur beam 33 .) From Fig. O one can infer that 



if the final temperature Tfi na j is identified with the hadron kinetic freeze-out temper- 
ature of, e.g., 125 MeV according to |2|, the corresponding initial temperature T initial 
would be 215 MeV. This statement, however, depends on the assumed temperature 
evolution. For instance, in case of a bag model equation of state a different weighting 
would occur and therefore different values of the initial and final temperatures would 
follow. Nevertheless, the merging of the initial and final temperatures, T initial and Tfl na i, 
at 170 MeV, displayed in Fig. [T^, suggests to use only T c s and N e g instead of more 
parameters. Indeed, the above time evolution equations and the values displayed in 
Fig. [12| reproduce the value of N eS = J dtV(t) as used in subsection 3.1. 



5 Summary 

An attempt is reported to explain nearly the whole set of the dilepton data in recent 
CERN-SPS experiments with heavy-ion beams. We assume that the dilepton emis- 
sivity is the same as that in qq annihilation. This hypothesis is made to work in the 
low mass region and in the intermediate mass region possibly for different reasons. 
Whether several channels add up to smooth out a spectrum with intrinsic structure 
( "kinematical saturation" ) or whether those structure are broadened and washed out is 
not addressed here. Note that while we find a good overall reproduction of the shapes 
of the experimental spectra, only the lead beam data can be explained with a unique 
and reasonable normalization. 

As main result we emphasize that the average temperature of T e fj = 160 • • ■ 170 MeV 
used in this study is in good agreement with the temperature deduced from measured 
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hadron ratios JT|. Since T e g is an average, one can conclude that the achieved maximum 
temperatures are above this value and, therefore, in a region where deconfinement 



is expected according to lattice QCD calculations f4"I |, which presently advocate a 
deconfinement temperature of T c = 170 MeV or somewhat larger. 

In the present work we restrict ourselves to central collisions (except for the NA38 
data) and compare several experiments with e + e~ and channels. On the basis 

of this study an analysis of the E± dependence of the combined NA38 and NA50 data 
seems feasible. This is of interest since, due to the impact parameter dependence, some 
interpolation from lead beam data to sulfur beam data is desirable. 

Our study is not a substitute for a detailed dynamical analysis. It is rather meant 
as a baseline calculation, designed to point out the underlying and unifying features 
of the experimental data, and to eventually extract a simpler physical message. We 
encourage dynamical microscopic prescriptions to attempt the global study performed 
here. 

With respect to the recent heavy-ion experiments at CERN-SPS with lower beam 
energies (40 AGeV) an interesting question is whether the featureless qq spectrum is 
compatible with the upcoming data. 
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A Background contributions 

Since we compose the dilepton spectra of incoherently added thermal and background 
contributions one has carefully to check the reliability of the background estimates. 
We employ here the event generator PYTHIA |26[] (version 6.104) for pp collisions and 
scale the results to AB collisions. Particular care has to be taken with the K factors. 
We use the parton distribution function set MRS D-'. 
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A.l Drell-Yan 



In Fig. a data compilation of the Drell-Yan cross sections J27| as a function of the 



scaling variable M/yfs and our PYTHIA simulations are displayed. The appropriate 
K factor is /Cdy — 1-2. Notice, however, that these data samples are for different beam 
energies and display a slight breaking of the anticipated scaling. Instead, independent 
fits of the invariant mass distributions of the data from |27|, ^ deliver an averaged K 
factor of 1.5. 

The intrinsic transverse parton momentum distribution can be fixed by the dilepton 
Qx spectrum in the Drell-Yan region M > 4.2 GeV. A comparison with the NA38 data 
[ p5| 1 delivers a value of (k ± ) = (0.8 GeV) 2 (see Fig. |TJ]). We use in all calculations (also 
for charm) this value. The K factor is less affected by variations of (k±). 

A. 2 Open charm 

We use here PYTHIA with charm mass parameter m c = 1.5 GeV and default fragmen- 
tation "hybrid". There are two different sources of a determination of the open charm 
K factor. Either one uses the hadronic channels, i.e. cross sections of identified D° 
and D + mesons from Fermi lab experiments, which are compiled in P5| , or the NA50 



dilepton data |30| for the reaction p(450 GeV) + W. Fig. [15| shows the comparison 
with hadron channels, which deliver K factors of JCpo = 5.4 and Kd+ — 6.8 resulting 
in an averaged open charm K factor of K, c = 5.7. Obviously these data do not con- 
strain the K factor very reliably. A sharper constraint is given by the dilepton channel 
which confirms this value as seen in Fig. [16]. The uncomfortably large K factor points 
to higher order processes which, however, could change the spectral shapes. In this 
respect the envisaged experiments by the NA60 collaboration [[!]]] to identify explicitly 
open charm are very important for this dilepton source. 
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Figure 1: Comparison of our model with the preliminary CERES data [pOR . Upper 
middle panel: invariant mass distribution, lower four panels: Q± spectra for various 
invariant mass bins (dashed lines: hadronic cocktail, solid curves: thermal yield, gray 
curves: sum). 
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Figure 2: Comparison of our calculations with NA50 data [^]. Left panel: invariant 
mass distribution, right panel: Q± spectrum for the invariant mass bin M = 2.1 • • -2.7 
GeV (solid curves: thermal contribution, dot-dashed curves: open charm, dashed 
curves: Drell-Yan, thin lines: parametrizations of the J/ip and ip' contributions ac- 
cording to [p5|] , gray curves: sum of these contributions). 




Figure 3: Coverage of the rapidity and transverse mass Mj_ of the various dilepton 
experiments, (a) left panel: lead beam, (b) right panel: sulfur beam. The arrows 
indicate the pp center-of-mass rapidity. 
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Figure 4: A Gaussian fit to the NA35 data of negatively charged hadrons [32]. The 
arrow indicates the center of the Gaussian which has a width of o^- =1.1. 
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Figure 5: Comparison of our model calculations with the CERES data JHJ with 
separately adjusted normalization factor N e g = 11.2 x 10 4 fm 4 (left panel) or with 
normalization N e g = 5.3 x 10 4 fm 4 adjusted to HELIOS/3 data [|37| (right panel). 
Dashed curves: hadronic cocktail from [Q, solid curves: thermal yield, gray curves: 
sum. 



17 



10 

0.5 1 1.5 2 2.5 3 3.5 
P T [GeV] 



Figure 6: A comparison of the thermal photon spectrum (without transverse matter 
flow) with the experimental upper bounds |34| when adjusting the source strength to 
the HELIOS/3 data " 




Figure 7: Comparison of our model calculations with the NA38 data |35] with sep- 
arately adjusted normalization factor N e g = 21.3 x 10 fm (left panel) and with 
normalization N e g = 5.3 x 10 4 fm 4 adjusted to HELIOS/3 data P7| (right panel). 
Meaning of the curves: thermal yield, open charm contribution and Drell-Yan (from 
left to right at larger M); gray curves: sum of all contributions; 
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Figure 8: Comparison of our model calculations with the transverse momentum 
spectrum of dileptons from NA38 |35j in the intermediate-mass region M — 2.1 • ■ -2.7 
GeV. Normalizations as in Fig. [7[ 
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Figure 9: Comparison of our model calculations with the normalized HELIOS/3 
difference data [S(200 A-GeV) + W minus p(200 A-GeV) + W] in various rapidity- 
bins with a unique normalization factor iV e fj = 5.3 x 10 4 fm 4 . The black dots 



indicate data groups for which also Q± spectra are available (cf. Fig. |10l) . 
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Figure 10: Comparison of our model calculations with the HELIOS/3 difference data 
[37]] for the M±_ distributions with a unique normalization factor N e g = 2 x 5.3 x 10 4 fm 4 
within the rapidity bin Y[ a b = 3.0 ■ ■ • 4.4 for the invariant mass intervals M = 0.2 ■ ■ ■ 0.6 
GeV (left panel) and M = 1.35- ■ -2.5 GeV (right panel). The additional factor 2 in 
the normalization is from a mismatch of these data compared with the ones in Fig. || 




2 3 
M [GeV] 



10" 



10" 



>" .,-.-11 

"o 10 



§ 10" 12 



10" 



10" 




12 3 4 

M T [GeV] 



Figure 11: Sensitivity of the dilepton spectra at midrapidity against transverse flow 
within the blast wave model. Left panel: dN/dM dY , right panel: dN/dM±dY . Solid 
curves: without transverse flow, dashed curves: with transverse flow v± = 0.6 (left 
panel) or 0.4 (right panel). The black (gray) curves are for a 2 (3) dimensional expan- 
sion (cf. |11[] for further details). 
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Figure 12: The change of the initial and final temperatures as a function of the 
normalization factor N. 
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Figure 13: A comparison of the Drell-Yan cross section delivered by PYTHIA (with 



K factor as described in text) with the data [27 . 
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Figure 14: Unnormalized transverse momentum distribution of dileptons in the Drell- 



Yan region, M > 4.2 GeV, in the NA38 experiment p5 |. The solid curve is our result 
from PYTHIA simulations with (k 2 ± ) = (0.8 GeV) 2 . A value of (k 2 ± ) = (1 GeV) 2 , 
as implemented in the recent PYTHIA version 6.143, results in an apparently worse 
description of the data. 
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Figure 15: A comparison of the open charm cross section (data compilation from |29|| ) 
with our PYTHIA calculations (with K factors as described in text). 
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Figure 16: A comparison of our PYTHIA calculations with the dilepton spectrum 



in p(450 GeV) + W reactions PD| . Dashed curve: Drell-Yan, dot-dashed curve: open 
charm contribution, gray curve: sum of all contributions. Included are parametrizations 
of the J I if) and if)' contributions according to pi] . 
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